DAMTP-200 7-4 MIFP- 07-02 
|hep-th/0701190| 

January 2007 



Gravitational Solitons and the Squashed Seven-Sphere 

P. Bizori 1 , T. Chmaj 2 ' 3 , G.W. Gibbons 4 and C.N. Pope 5 

1 M. Smoluchowski Institute of Physics, Jagiellonian University, Krakow, Poland 
2 H. Niewodniczanski Institute of Nuclear Physics, Polish Academy of Sciences, 

Krakow, Poland 
3 Cracow University of Technology, Krakow, Poland 
A D.A.M.T.P., Cambridge University, Wilberforce Road, Cambridge CBS OWA, UK 
5 George P. & Cynthia W. Mitchell Institute for Fundamental Physics 
Texas A&M University, College Station, TX 77843-4242, USA 



ABSTRACT 

We discuss some aspects of higher-dimensional gravitational solitons and kinks, includ- 
ing in particular their stability. We illustrate our discussion with the examples of (non-BPS) 
higher-dimensional Taub-NUT solutions as the spatial metrics in (6 + 1) and (8 + 1) di- 
mensions. We find them to be stable against small but non-infinitesimal disturbances, but 
unstable against large ones, which can lead to black-hole formation. In (8 + 1) dimen- 
sions we find a continuous non-BPS family of asymptotically-conical solitons connecting 
a previously- known kink metric with the supersymmetric A$ solution which has Spin (7) 
holonomy. All the solitonic spacetimes we consider are topologcally, but not geometrically, 
trivial. In an appendix we use the techniques developed in the paper to establish the linear 
stability of five-dimensional Myers-Perry black holes with equal angular momenta against 
cohomogeneity-2 perturbations. 
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1 Introduction 

1.1 Gravitational solitons and kinks 

By analogy with other areas of physics, a Gravitational Soliton in n spacetime dimensions 
may be defined to be an everywhere complete non-singular globally stationary Lorentzian 
spacetime M, satisfying the vacuum Einstein equations [ljo Thus a gravitational soliton 
has by definition neither an event horizon nor an ergo-sphere and should therefore be dis- 
tinguished from a stationary or static black hole, which is only required to be non-singular 
outside a a regular event horizon. For conventional solitons in flat space, one usually adds as 
a requirement that the solution not only be non-singular, but also have finite total energy. 
Furthermore, this energy is determined by the conserved charges the soliton may carry and 
also by quantities specifying asymptotic boundary conditions. Thus, for example, a mag- 
netic monopole in Yang-Mills theory has a mass determined by its magnetic charge and by 
the vacuum expectation value of the Higgs field at infinity. 

There are no gravitational solitons without horizons in four spacetime dimensions but 
in higher dimensons such objects do exist: perhaps the best known example being the 
Kaluza-Klein monopole in five dimensions [21 [7] . This has a metric of the form 

dsl = -dt 2 + ds\ , (1) 

where ds\ is the self-dual Taub-NUT gravitational instanton metric. This four-dimensional 
metric is asymptotically locally flat, with a circle direction that stabilises to a constant 
length at large distance. The Kaluza-Klein monopole, upon reduction along this circle to 
(3 + 1) spacetime dimensions, has a finite ADM mass given by the length of the circle at 
infinity. 

In this paper we shall be concerned only with the case of vanishing cosmological constant although many 
of the ideas go through in the case that the cosmological constant is negative 
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1.2 Vacuum interpolation 

In common with many flat-space solitons, the Kaluza-Klein monopole may be thought 
of as spatially interpolating between two inequivalent "vacua" or "ground states" of the 
theory, namely the flat five-dimensional Minkowski spacetime E 4,1 near the origin, and the 
compactified Kaluza-Klein ground state S 1 x E 3 ' 1 near infinity. If a solitonic solution of the 
Einstein equations exhibits vacuum interpolation, it seems reasonable to refer to it also as 
a Gravitational Kink^ 

Again by analogy with other areas of physics, one does not expect to find an Asymptot- 
ically Flat gravitational soliton, i.e. one which outside a compact spatial set or world tube 
tends to the flat metric on n-dimensional Minkowski spacetime E n_1,1 . This is because it 
would interpolate between two copies of the same vacuum. Indeed, as we shall see shortly, 
this intuitive expectation is borne out by a No-Go theorem. Thus with respect to the flat 
Minkowski vacuum E n ' , one cannot think of a gravitational soliton as having finite en- 
ergy with respect to the flat Minkowski vacuum. Nevertheless, with respect to the ground 
state near infinity it may well have finite energy. 

Vacuum interpolation also occurs in the case of extreme black holes or extreme black 
p-branes. Again, this is between different kinds of ground states; typically between a flat 
vacuum at infinity, and a compactified AdS p +2 x K, where K is an (n — p — 2)-dimensional 
compact manifold [5]. 

1.3 Classical stability 

The definition given above does not specify whether a gravitational soliton should be stable. 
That is deliberate, because although to be important as a potentially long-lived classical 
"lump," to use Coleman's phrase [6], a gravitational soliton should at least be be classically 
stable against small but non-infinitesimal disturbances (i.e. not merely linearized fluctu- 
ations). It is not reasonable, however, to demand classical stability, in any gravitational 
theory, against all possible large disturbances, since nothing forbids gravitational collapse 
to a black hole with the same asymptotics. After all, even Minkowski spacetime is unstable 
against the possibility of a concentrated region of gravitational waves undergoing gravi- 
tational collapse to a black hole. Indeed in a recent numerical study of the dynamics of 
Kaluza-Klein monopoles [TO], precisely such a collapse was seen to occur. The possibility 
of black hole formation, and the resulting spacetime singularities, invalidate the type of of 
topological stability criteria derived from cobordism theory that were discussed in 

Both Minkowski spacetime and the Kaluza-Klein monopole are supersymmetric, or BPS. 
Thus these two examples clearly demonstrate the fallacy of the widespread belief that to 
prove stability in general relativity it suffices to establish that the spacetime admits a Killing 
spinorH 

The singularity theorems of classical general relativity show that these black holes con- 
tain spacetime singularities, which are a clear indication that the classical theory is incom- 

2 Note that the use of the word kink here should be distinguished from the notion of the "kink number" 
of a Lorentzian metric with respect to some hypersurface, introduced by Finkelstein and Misner [3], and 
elaborated upon in [3]. 

3 Quite apart from the non-linear instabilities involving black-hole formation, further instabilities of 
Minkowski spacetime, or indeed any asymptotically flat spacetime, can arise unless one restricts attention 
to perturbations or deformations that decay at large distances. The example of Kasner spacetime clearly 
shows that flat space is unstable to the formation of all-encompassing naked singularities in finite time, if 
one allows perturbations that do not decay near infinity. 
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plete. If these singularities are hidden inside event horizons, i.e. if cosmic censorship holds, 
they may not be an obstacle to studying the evolution of the exteriors of black holes. How- 
ever, such black holes do not have a fixed mass, and may grow by, for example, absorbing 
radiation. Thus in general, black holes cannot be thought of as solitons. An exception may 
arise if one considers so-called extreme black holes, in which the mass may be determined 
entirely in terms of conserved charges [12]. We shall not discuss this type of "solitonic" 
black hole further in this paper. 

1.4 Quantum stability 

The stability considerations described above were purely classical. Quantum mechanically, 
black holes are unstable against thermal Hawking radiation. Thus in the case of the col- 
lapsed Kaluza-Klein monopole, it seems very likely that the ultimate quantum-mechanical 
state is the monopole itself, since the magnetic charge cannot be radiated away. 

If, as is currently widely believed, the evaporation of neutral black holes leads to their 
complete disappearance, it would seem that Hawking evaporation is essential in order to 
solve the problem of classical singularities. 

A frequently used criterion for the stability of a particle in quantum mechanics is that 
it has the least mass of any state carrying the same conserved charges. This criterion, 
often associated with BPS bounds, is often used to argue that various solitons, or indeed 
ground states, are stable. In the case of spacetimes, what is often in one's mind is quantum 
tunnelling. While it is certainly true that a BPS condition, or the existence of a Killing 
spinor, may mean that tunnelling is impossible, it does not rule out the sort of classical 
instabilities we discussed previously. 

1.5 Ultra-Staticity 

The main subject of the present paper is the case of gravitational solitons in nine spacetime 
dimensions. This dimension is large enough to admit a rather richer structure of solitonic 
solutions than can be obtained in lower dimensions. Before describing our new results 
however, it may be useful to continue the general discussion, making it a little more precise. 
In particular we wish to establish the general result that that a gravitational soliton, as we 
have defined it, must be ultra-static, i.e. the (unwarped) product M = E x S of time with 
a complete non-singular Ricci flat spatial Riemannian (n — l)-manifold E. 

The assumption of global stationarity means that the spacetime is an R-bundle over £ 
the space of orbits of time translations. The metric may thus be cast in the form 



where i = 1,2, ... ,n — 1, and the everywhere non-vanishing strictly positive function V 
and the Sagnac R-connection ujidx 1 are independent of time. The curvature of the Sagnac 
connection is given by 




(2) 



Fij — diiOj — djLOi . 
We begin by noting that the vacuum Einstein equations imply that 



(3) 




(4) 



Multiplication by ujj and integration by parts gives 
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If the boundary term at infinity vanishes, then we conclude that 

F lj = 0. (6) 

The Einstein equations then imply 

VV = 0. (7) 

Thus if V is bounded at infinity, a standard argument based on the maximum principle 
shows that 

V = constant (8) 
The remaining Einstein equation then implies that the spatial metric is Ricci flat 

Rij = 0. (9) 

If we assume that £ is simply connected, or pass to a finite covering space if it is not, we 
may set Ui = 0, V = 1 and the metric is ultra-static 

ds 2 = -dt 2 + g lj dx i dx j . (10) 

It follows from (|10j) that the question of the existence of gravitational solitons reduces 
completely to that of the existence of complete Ricci-flat spatial manifolds {£, <7y}. It is 
known [9] that there are no non-trivial Asymptotically Euclidean^ metrics, and hence no 
asymptotically flat gravitational solitons but there are plenty of metrics which are Asymp- 
totically Locally Euclidean^ as well as metrics with much more complicated asymptotics. 



1.6 BPS solitons 

Among the various possibilities for the Ricci-flat spatial metric, of particular interest are 
those admitting a covariantly constant spinor. They have reduced holomony, and if n < 
12, the spinor field is a Killing spinor of a supergravity theory and the solitons are thus 
supersymmetric. The existence of the Killing spinor allows one to relate the spectrum of the 
Lichnerowicz operator acting on symmetric traceless second-rank tensors to the spectrum 
of other differential operators. In this way, one may establish the linearized stability of 
solitons with special holonomy. For example, a metric with Spin (7) holonomy admits a 
covariantly constant self-dual 4-form. This 4-form may be used |20| to establish a 1-1 
correspondence between the spectrum of the Lichnerowicz operator and the spectrum of 
the Hodge de-Rham operator acting on anti-self-dual 4-forms. Since the spectrum of the 
latter is manifestly non-negative, it follows that the Lichnerowicz operator has no modes 
of negative eigenvalue, and hence that Spin(7) solitons are stable at the linearised level. 
However, as we discussed earlier, they will not be stable against deformations sufficiently 
large that collapse to black holes takes place. Similar remarks about linearised stability 
apply to the other cases of special holonomy, which circ cis follows: 

• Ricci-flat Kahler or Calabi-Yau with holonomy SU(k) C SO(2k) and thus n = 2k + l, 

• Hyper-Kahler with holonomy Sp(k) C 50(4fc), and thus n = 4k + 1, 

• Holonomy £?2 C SO (7), and thus n = 8, 

4 i.e. which tend to the flat metric on (n — l)-dimensional Euclidean space E"" 1 outside a compact set 
5 i.e. which tend to the flat metric on E n_1 /r , F C 0(n — 1) outside a compact set 
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• Holonomy Spin(7) C SO(8), and thus n = 9. 

Explicit complete non-singular metrics are known in all cases. The easiest examples to 
construct are cohomogenity one and are Asymptotically Conical(AC); they tend to Ricci 
flat cones over an Einstein manifolds which are: 

• Holonomy SU(k) C SO(2k): Einstein-Sasaki, 

• Holonomy Sp(k) C SO (4k): Einstein- Tri-Sasaski, 

• Holonomy G 2 C SO (7) : Weak SU(3), 

• Holonomy Spin(7) C 50(8): Weak G 2 . 

Other explicit cohomogenity one examples have been found which are Asymptotically 
Locally Conical (ALC). In this case a circle subgroup of the isometry group has orbits 
which tend to constant length at infinity. The ur-example is the Taub-NUT metric, i.e. the 
Kaluza-Klein monopole. 

1.7 Cohomogeneity One and Cohmogeneity Two 

In this paper we shall restrict attention to complete Ricci flat (n — l)-dimensional positive 
definite metrics which are of cohomogeneity one, that is whose isometry group G has prin- 
cipal orbits which are (n — 2) dimensional. Such solutions give rise to static solitions on 
n-dimensional Lorentzian spacetime obtained by taking the metric product with time. The 
isometry group of the n-dimensional static spacetime is therefore the product IxG. 

We then construct a consistent time-dependent ansatz for the n-dimensional spacetime, 
which is is invariant under just the action of the orginal group G, and which agrees with 
the static soliton solution in the special case that there is no time dependence. The general 
time-dependent Lorentzian spacetime is thus of cohomogeneity two. The word "consistent" 
means that every solution of the resulting system of (1 + l)-dimensional equations gives a 
solution of the n-dimensional vacuum Einstein equations. 

The reason for restricting to cohomogeneity one and cohomogeneity two is not only for 
simplicity but because it allows us to exploit the considerable body of existing information 
in the literature on cohomogeneity-one Ricci- flat metrics. 

We shall also mainly concentrate on the case when the spatial manifold is topologically, 
but not geometrically R n . 

2 Higher-Dimensional Time-Dependent Taub-NUT Solitons 

One may consider a variety of higher-dimensional static metrics of the form ds 2 = —dt 2 +ds 2 , 
where ds 2 is a Ricci-flat spatial soliton metric. Then, following the same strategy as in 
[TU] , these metrics may be used to provide initial data for the time-dependent vacuum 
Einstein equations. Following the discussion in [TO], we shall consider the higher-dimensional 
analogues [131 [14] of the four-dimensional self-dual Taub-NUT metric. It should be noted, 
however, that these higher-dimensional analogues are non-supersymmetric, in the sense that 
unlike the four-dimensional Taub-NUT case, there is no Killing spinor. We shall discuss the 
examples of the Taub-NUT metrics in six and eight dimensions below, after first presenting 
a general time-dependent ansatz. 
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2.1 The time-dependent ansatz 

A suitably general time-dependent ansatz for our purposes is 

d f = _ Ae - 25 dt 2 + A- 1 dr 2 + r 2 [e 3B dY? m + e~ 6mB a 2 } , (11) 

where dE^ is the metric on an Einstein-Kahler space of real dimension 2m, normalized so 
that it satisfies R a b = 2(m + l)g a b- The 1-form a is given by a = dtp + B, where B is a 
potential on the Einstein-Kahler base space, such that dB = 2J, where J is the Kahler 
form. The total dimension of the spacetime is D = 2m + 3. 

Special cases of (fTT]) include, if e 2S = A and the metric is assumed to be independent 
of time, the Taub-NUT solitons. For these solutions, the Einstein-Kahler base spaces are 
taken to be CP m . Other important special cases are the higher-dimensional Kerr and Kerr- 
AdS metrics. In the case of (2m + 3) dimensions, with all rotation parameters set equal, 
these have cohomogeneity one, and their stability can be analysed by a small extension of 
our procedure in which the Kaluza-Klein vector is retained also in the reduction. This is 
discussed in the Appendix. 

The time-dependent Einstein equations Rmn = 2(m+l)A i^vfiV for the metric (fTT|) break 
up into momentum and Hamiltonian constraint equations 

, = _ 2mA 3B _ 6{m+l)B] _ l mT ,^ A -x b 2 + AB' 2 ) 

r (2m + l)r 1 V ' J 2 v ' 

„-4m-3 e 6mfl -2(m+l)Ar, 



2(2m+ 1) 

A = —9mrABB' , (12) 



a slicing condition 



and a wave equation 



S ' = -l mr (e 25 A' 2 B 2 + B' 2 ) , (13) 



e S A -l r 2m+l B \ ■ _ ( e S Ar 2 m +l B ,V + 4(m+l) s 2w _! ?jB _ 6{m+1)B) 
J V J 3(2m + l) V ; 

+ — — r ~ 2m ~ 3 e~ 5 e 6mB = . (14) 

2(2m + 1) y ' 

If we define a mass functional Mir, t), by writing A = 1 — M(r, t)/r 2m — Xr 2 , then the 
Hamiltonian constraint becomes 



,2 



M' = -mr 2m+1 (e ^4^ 1 i? 2 + AB ) + 



— 2m— 3 g6m_B 



2 v ' 2(2m + l) 

+ _^_ r 2m-l [(2m + 1} + e -6(m+l)B _ 2(m + 1)e -3B] (15) 

Note that the right-hand side is manifestly positive. 

In what follows, we shall specialise to the case of the vacuum Einstein equations, by 
setting A = 0. 
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2.2 The (6 + l)-dimensional Taub-NUT soliton 

The spatial metric of the time- independent (6 + l)-dimensional Taub-NUT soliton is given 
by 

2 (p + £) 2 dp 2 2 £ 2 (p-£)(p + 3£) , , 2 2 
= 2(p -*)(, + 3*) + ^TIP +(P -°^ 2 ' (16) 

with p > £, where the notation for a and cE 2 here is the same as in (jlip . One can see that 
the metric near p = I approaches the origin of hyperspherical coordinates in R 6 , by defining 
a new radial coordinate y = y/2£(p — £). At large p, the metric approaches R 5 times a circle 
of asymptotic radius y2£. The manifold on which (fT6j) is defined has the topology of R 6 . 

Comparing with the ansatz (jlip . with m = 2, we see that the six-dimensional Taub-NUT 
metric gives initial data with 

15Bo _ 2£ 2 (p + 3£) 2(p + 3£)(p-£) ( dry 2Sq 

where 

r w = 2£ 2 (p-£f(p + £) 2 (p + 3£). (18) 

2.3 The (8 + l)-dimensional Taub-NUT soliton 

The spatial metric of the time- independent (8 + l)-dimensional Taub-NUT soliton is given 
by 

, Hp + £f dp 2 8£ 2 (p - £)(p 2 + Up + 5l 2 ) 2 2 2 2 
dS8 = 8(p-£)( /) 2 + 4£p + 5£ 2 ) + +(P " } 3 ' (19) 

with p> I. It is defined on the manifold R 8 . 

Comparing with (jlip with m = 3, we see that the eight-dimensional Taub-NUT metric 
gives initial data with 

^ 21Bo _ 8£ 2 (p 2 + A£p + 5£ 2 ) _ 8 (p-£)(p2+4£p + 5£ 2 ) ( dr^ 2Sq _ 

- ^pTW ' W+W \~d~p) ' " 

(20) 



where 



r 14 = | £ 2 (p + £f{p - £) 7 {p 2 + 4£p + U 2 ) . (21) 
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As we shall see below this solution is a special case of a more general ansatz. 



3 Nine-Dimensional Squashed Seven-Sphere Solitons 

Our previous ansatz (jlip was based on the Hopf fibring of S 2m+1 by U{\) Hopf fibres 
over a CP m base manifold^. In the case where m = 2p, one can instead consider S 4p+3 
regarded as an 577(2) bundle over HP P . The simplest case is for p = 1, with S 7 regarded 
as an SU(2) bundle over S 4 . In this section, we shall consider a time-dependent ansatz for 
an (8 + l)-dimensional time-dependent metric where the spatial 8- metric has surfaces at 
constant radius that have S 7 topology, fibred by S 3 . Two deformation parameters will be 
included, one parameterising the volume of the S 3 fibres, and the other parameterising a 
homogeneous squashing of the S 3 fibres themselves. 

6 Or fibrations over more general Einstein-Kahler base manifolds. 
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We begin with some group-theoretic preliminaries, by considering left-invariant 1-forms 
Lab for 50(5), These obey Lab = —Lb A and 



cILab = Lac A Lqb 



(22) 



We take the 50(5) indices to range over < A < 4, and split them as A = (a, 4), with 
< a < 3. The 50(4) 1-forms L a b are then expressed in an SU(2)l x SU(2)r basis with 
generators 

Ri = \{L$i + \eijkLjk) , Li = |(Loi — ^tijk-Ljk) > (23) 
where 1 < i < 3. The 5 7 = 5O(5)/50(3) coset is then spanned by the 1-forms 



Pa — \L a ^ ■ 



(24) 



(Note a rescaling of P a , relative to [16]. This is done for convenience, to avoid factors of 2 
later.) The algebra of the 1-forms is easily seen to be 



dP a 
dLi 



—tijkRj A Rk — J l ah P a A Pb- 

Jab Ri A A + Jab L% A P& , 

e ijkLj A Lfc — P a A Pft , 



(25) 



where we have defined antisymmetric self-dual and anti-self-dual 'tHooft tensors J l ab and 

J l ab by 



The metric on the unit round 5 7 is given by 



J 



jk ~ t ijk 



J 



jk 



^ijk ■ 



dti 2 . 



Rf + p; . 



(26) 



(27) 



Our general time-dependent ansatz is 

dsl = -A e~ 2S dt 2 + A- 1 dr 2 + r 2 {e" 4B [e 2C (R 2 + R 2 ) + e" 4C R 2 } + e 3B P 2 } . (28) 

Straightforward calculations show that the Ricci-flatness of d§g implies the following equa- 
tions. First, we have the Hamiltonian and momentum constraints 

A' = 



+ 



6 . 6r , 

-A A 

r 7 

2 r 



-(B' 2 + e 25 A' 2 B 2 ) + (C 2 + e 2S A~ 2 C 2 ) 



7r 



e 4B-8C , 4e 4B-2C _ 4g - 10B+2C _ 2g - 10S-4C 



+ 24e 



-3B 



.4 



-6rA(P^ / + -00 / ). 



In addition, there is the slicing constraint 
5' = -3r 



{B > 2 + e 2<5 ^-2 £2) + 2 (c ,,2 + e 25 A -2 £ 2) 



(29) 

(30) 
(31) 
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Finally, we have the dynamical equations for the two squashing modes, which give 
eV A' 1 r 7 B) - e s {e~ 5 A r 7 B')' 

(_ e 4B-8C + 4e 4B-2C + 10e -10B + 2C + 5e ~10i?~4C_ 18e -3B ) = ^ ^ 



e V A" 1 r 7 C) - e\e~ & A r 7 C')' 

+ ^1 (_ e ^-8C + e 4B-2C + e -10B+2C _ e -10B-4 C) = q (33) 

3 

It can be straightforwardly verified that the constraints are indeed consistent with the 
dynamical equations of motion. 

As a check, it can be verified that if we set C = ^B, for which the ansatz (f28l) reduces 
to the special case of setting m = 3 and A = in (fTTj) . i.e. describing a squashing of S 7 
viewed as a U(l) bundle over CP 3 , we indeed obtain the same equations as those given in 
section 12.11 Another check is instead to set C = 0, in which case the system reduces to the 
one discussed in [T7], where S 7 is viewed as a round S 3 bundle over S 4 . 

3.1 Static solutions 

In this section we consider regular static asymptotically (locally) conical solutions of the 
system (29-33). Note that in agreement with Section 1.5 all static solutions are ultrastatic, 
i.e., Ae~ 2S = 1, thus equations (29-33) reduce to the following system of ordinary differential 
equations 

e 5 (e s r 7 B'y + — (- e ^-8C +4e 4B-2C + 10e -ios+2C + 5e -iOB-4C_ 18e -3B ) = 0> (34) 

e s (e s r 7 Cy - — {-e iB ^ c + e 4B ~ 2C + e - 10B+2C - e~ WB - iC ) = . (35) 
3 

5' = -3r {B' 2 + ^C' 2 ) . (36) 
Regularity at the origin implies the following behavior for r — > 

B ~ br 2 , C ~ cr 2 , e 5 = 1 - (3b 2 + ^c 2 )r 4 , (37) 

where b and c are free parameters. Using scaling symmetry, without loss of generality, we 
can set 6 = 1. Then, (|37p gives rise to a unique one-parameter family of local solutions 
parametrized by c. Numerical analysis shows that for any c in the interval < c < 7/2, 
these local solutions can be continued to infinity and thus give the desired global solutions 
(see Fig. 1). For the values of c outside this interval the solutions become singular for a finite 
r. The asymptotic behavior of global solutions near infinity depends on c: for < c < 7/2 
the squashing modes grow logarithmically and 6 goes to minus infinity, while for c = both 
B and 5 have finite limits. Thus unless c = 0, the space sections are asymptotically locally 
conical, ALC, but in the limiting c = case they become asympototically conical (AC). For 
two values of c the solutions are known is closed form: these are the Taub-NUT solution 
(|19p which corresponds to c = 7/2 and the so called Ag solution [151 CE] which corresponds 
to c = 2. 

Below we discuss in detail the structure of static solutions and their stability properties. 
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Figure 1: Plot of C/B for static solutions for several values of the shooting parameter c. For 
< c < 7/2 we have C - (7/2)5 - (1/2) In 2 for r -> oo. 



3.2 The Spin(7) background A 8 (c = 2) 

The As solution is an ultra-static nine-dimensional vacuum solution whose space sections 
are asymptotically conical (AC). 

ds 2 9 = -dt 2 + dsj , (38) 

where ds 2 , is a Ricci-flat metric with Spin(7) holonomy. A simple metric of this kind, which 
extends smoothly onto a manifold of R 8 topology, was obtained in [151 116] . where it was 
denoted by As. In the normalisation we are using here, it is given by 



2 (p + tyd P 2 2 2 At\p + u){ P -i) 2 2 

ds 8 = ( p + 3£)(p_^ + ^ + 3 ^^~ £ )( i? l +jR 2) + (p + l)2 ^3 + 2^ ) P a • (39) 

The radial coordinate p lies in the range £ < p < oo. Near p = £ we may define a new radial 
coordinate y = 2\J (p — £)£, in terms of which the metric approaches 

ds 2 ~dy 2 + y 2 (R 2 + P%) (40) 

at small y. At large distance, p —* oo, the metric approaches approaches R 7 times a circle 
of asymptotic radius 21. The situation is therefore closely analogous to that of the self-dual 
Taub-NUT metric in four dimensions. 

Expressed in terms of the ansatz (|28p . the As solution has the form 



„21B 



2{p + tf , c _ 21 



e = "~-t7tt , e 



^ 2 (/9 + 3£) 3 ' + 
(p + 3l)(p-l) /drx2 
(p + £) 2 \dp) ' 



where 



e 25 = A, (41) 

\p-r<-)- ^upy 
r 14 = 64£ 2 {p + 3£) 3 (p - £) 7 (p + £f . (42) 
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3.3 The continuous family of solutions (0 < c < 7/2) 

We define the new independent variable r by re s d/dr = d/dr and let x = B and y = y| C. 
Then, assuming staticity, equations (f34"|) and (f35|) take the form 

d 2 x s dx dV d 2 y sdy dV , A , 

— + 6e 5 — + — = 0, -4 + 6e <5 / + — = 0, 43 

where 

V (x, y) = ^- ( - e 4 ^" 4 ^ + 4 e ^-^ _ ^-iox+v^, _ 2e -i0x- 2v /i% + 24e -3^ ^ (44) 
and the slicing constraint (|36p becomes 

^ J - 3 ((l) 2 + (f) 2 )- < 45 » 

The boundary conditions (|37p translate to the following asymptotic behavior for r — ► — oo 

x~6e 2T , y~y|ce 2T , e 5 ~ 1 - 3(6 2 + ^c 2 )e 4r . (46) 

As above we set b = 1, hence we have a one-parameter family of local solutions parametrized 
by c. 

It is useful to interpret the above system in terms of the mechanical analogy of a sticky 
ball rolling on the surface z = V(x,y). Due to the friction the energy of the ball, 



E = ~ 
2 



\drJ \drJ 



+ V(x,y), (47) 



dE ~j((dxy fdyy 



dr ^ W + B^°' 



decreases in time 



Note that the combinations of equations (j45J) and (|48p together with the boundary condi- 
tions (|46p yield the constraint E = e 2S which can be used to eliminate e* 5 from equations 

Assuming that x and y are both positive and large we can solve the equations (|4"3"]) -(|45 p 
asymptotically to get the following two possibilities 

/7 ,. .x Ff hi 2 

(i) y~\j 2 x or w y~^2 x ~ ~7u ' ^ ' 



The case (i) is exceptional and corresponds to the ball rolling down the ridge, while the 
case (ii) is generic and it corresponds to the motion down the valley (on the left side of the 
valley there is the ridge which separates it from a cliff and on the right side there is a steep 
ascent). In the second case the asymptotic behavior of solutions for r — > oo is 

x~-mr In , e ~— . (50) 

3 3 V 10 / r y ' 
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Figure 2: The potential z — V(x,y) and three trajectories which start at the peak (1,0,0) with 
different slopes (that is, different values of the parameter c). The black curve (c — 7/2) represents 
the Taub-NUT solution, the dark blue curve (c < 7/2) represents a generic solution of the continuous 
family, and the red curve shows a singular trajectory with c > 7/2. 

3.4 Stability of static solutions 

The role of static solutions in dynamics depends on their stability properties. In this section 
we investigate the stability of static solutions presented above (0 < c < 7/2). 

3.4.1 Linear stability 

Following the standard procedure we seek solutions in the form 



where the index denotes a static solution and the index 1 denotes a perturbation. We 
substitute this expansion into equations (29)- (33) and linearize them. Integrating equation 





(51) 
(52) 
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(30) we obtain 
and from equation (31) we get 



A l = -§rA G {B' Q B 1 + ^C' C 1 ), (53) 



S[ = -6r(B' B[ + ^C' C[) . (54) 

Inserting equations ([53]) and ([54"J) into the linearized equations (32) and (33) and separating 
the time dependence B\{t,r) = ex.p(—i\t)v\(r), C\{t,r) = ex.p(—i\t)u\(r) we get the 
eigenvalue equation for the spectrum of small perturbations 

-l e 5 °(>r 7 S A y + KS A = A 2 S A , (55) 

where 



(56) 

and K is a 2 x 2 matrix determined by the static solution. The matrix K is very complicated 
but fortunately we do not need it in an explicit form. To demonstrate stability we exploit 
the existence of the zero mode 



C>(r)J ' (57) 

which is due to the scaling invariance of the problem. From the above heuristic analysis 
of the behavior of static solutions it is clear (although not proved rigorously) that the 
functions -E>o( r ) and Co(r) are monotone increasing, hence the components of the zero 
mode have no zeros. Thus, it follows from the Sturm-Liouville theory that there are no 
negative eigenvalues. 

We point out that in the case of the A$ manifold, the linearized stability of this solu- 
tion follows from the existence of a covariantly-constant spinor in the Spin(7) holonomy 
manifold As, as discussed in section 11.61 Specifically, the spectrum of the Lichnerowicz 
operator describing transverse traceless metric perturbations is identical to the spectrum 
of the Hodge-de Rham operator acting on anti-self-dual 4- forms [20J, and thus there can 
be no negative-eigenvalue modes and hence the solution is stable at the level of linearized 
perturbations. 



3.4.2 Nonlinear stability 

In order to verify numerically the nonlinear stability of static solutions numerically, we have 
expressed equations (29)- (33) in the first order form using the momentum variables 

p B = ^A^B, P c = e 8 A~ x C . (58) 

We have solved the resulting equation system using the free evolution scheme in which the 
function A(t, r) is updated from the momentum constraint (30). Integration in time is done 
by the modified predictor-corrector McCormack method on a uniform spatial grid. The 
slicing constraint (31) is solved with the fourth order Runge-Kutta method. The whole 
procedure is second order accurate in time and fourth order in space. The results shown 
below were produced for initial data of the form 

B(0,r) = B (r), C(0,r) = C (r), P B (0, r) = p ^ e"T, p c ( , r ) = 0, (59) 
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where the amplitude p was varied and the parameters R and s were kept fixed. We have 
found that for small perturbations, that is for small values of the control parameter p, the 
solution returns to equilibrium and the excess energy of the perturbation is radiated away 
to infinity, while for large perturbations a black hole forms. The behavior is qualitatively 
the same for all static solutions (independently of c) and we illustrate it in Figs. 3 and 4 in 
the case of the As solution. 
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Figure 3: Asymptotic stability of the solution Ag. For initial data (|59|) with a small amplitude 
(p = 0.1, R = 3, s — 1) we plot a series of snapshots of the function B(t, r) (where t is central proper 
time). The dashed line shows the unperturbed Ag solution. 



3.5 The Eight-Dimensional Kink Solution (c = 0) 

In this section, we shall show the existence of a complete Ricci-flat 8-metric, which is 
asymptotically conical (AC), which we shall call the kink. It may be considered as spatially 
interpolating between flat Euclidean 8-space near the origin, 

B = 0, C = 0, A = e 25 = 1 , (60) 

and the Ricci-flat cone over the Einstein-squashed 7-sphere at infinity, 

B = ^, C = 0, A = e 2S = 9-5- 10 / 7 . (61) 

It was demonstrated in [18], at the numerical level, that there exists a complete and non- 
singular metric that interpolates between these two constant solutions. It is defined on a 
manifold of R 8 topology. Here, we present a sketch of a proof of the existence of this kink 
solution. 

Repeating the steps from Section 3.3 and setting y = in equations (f4"3j) - (f4"5j) . we obtain 
the 3-dimensional autonomous system 

d 2 x sdx dV , , 

^ +& V + ife= ' (62) 
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Figure 4: Instability of the solution As for large perturbations. For initial data ([59]) with a large 
amplitude (p — 0.3, R = 3, s = 1) we plot a series of snapshots of the function A(t, r) (where t is 
central proper time). During the evolution A(t, r) drops to zero at r = ru ~ 1.773 which signals the 
formation of an apparent horizon there. Outside the horizon the solution relaxes to a static black 
hole. 



where 

V = i (8e- 3x + e 4x - 2e- 1Cte ) . (64) 

The boundary conditions (|37|) for r — * — oo simplify to 

x~be 2r , e 5 ~ 1 - 36 2 e 4T , (65) 

hence up to scaling given by b we have a unique regular local solution. 

As above, we can interpret this system in terms of the mechanical analogy of a ball 
rolling in the potential V(x) with a variable friction (see Fig. 5). The energy of the ball 

*-5(i) ,+ "M < 66 > 

As in the general case we have the constraint E = e 2S . Using this constraint we eliminate e s 
from equation (j62[) and get the autonomous 2-dimensional dynamical system. This system 
has two critical points which correspond to the constant solutions ([60]) and (f6Tj) : the saddle 
(x = 0,dx/dr = 0) and the stable node (x = ^,dx/dr = 0). The function e 2S serves as 
the Lyapunov function, thus it is evident that the orbit starting from the saddle (0, 0) along 



decreases in time because 



16 



1.08 
1.04 
1 

0.96 
0.92 
0.88 
0.84 

-0.1 0.1 0.2 0.3 0.4 

V 

Figure 5: In terms of the mechanical analogy the kink solution corresponds to the ball rolling down 
from the maximum of the potential V — 1 at x — to the minimum V = 9 • 5 -10 ^ 7 ~ 0.903 at 

t — 1" 5 
x — 7 . 

the unstable manifold will end up at the stable node 0) for r — > oo. The linearization 
around this critical point yields the following asymptotic behavior for r — > oo 

(%7)=^ a,t «'+^ 

where the eigenvalues are 

A 1 = -^5 2 / 7 , A 2 = -2-5 2 / 7 , 
5 

and the corresponding eigenvectors are 

The kink orbit approaches the node (^,0) along the slow eigendirection £i. We claim 
that the kink orbit stays for all times in the first quadrant of the (x,dx/dr) plane. To 
see this, consider an exceptional trajectory which approaches the node (^,0) along the 
fast eigendirection £2 (that is, c\ = and C2 < in equation ([55]) ). This trajectory, run 
backwards in time from the node, obviously cannot cross the line x = ^ and consequently 
it prevents the kink trajectory (which is trapped below the fast eigendirection trajectory) 
to do so. In terms of the mechanical analogue of a ball rolling in the potential V the above 
analysis demonstrates that the motion of the ball is overdamped and moreover the ball 
rolling down from the maximum of the potential at x = cannot overshoot the minimum 
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Figure 6: The orbit of the kink. The node is approached along the slow eigendirection £i. 
3.5.1 Analytic study of the kink solution 

In this Section we describe an attempt to find the kink solution in the closed form. Although 
this attempt was not successful, we believe it is worth presenting because it yields a deeper 
analytic insight into the structure of equations. In addition, it provides an alternative way 
of proving the existence of the kink. 
Consider the 8-metric 



dsg = dr 2 + a 2 (<Tj 



a 1 ) 2 + b 2 dn 2 , 



(71) 



and a 



where dQ, 2 is the metric on the unit 4-sphere, A 1 is the 1-instanton solution on the S 4 , 
are the left-invariant 1-forms of SU{2). The flat metric on M 8 corresponds to a = b = ^r, in 
which case the principal orbits are round 7-spheres for all values of r. In the kink solution, 
which is easily found by numerically solving the Ricci-flat equations for the metric ansatz 
(|7ip . the metric approaches the flat form at small r, whilst as r goes to infinity the principal 
orbits approach the squashed Einstein metric on S 7 , for which a 2 /b 2 = 1/5. In fact, as r 
approaches infinity the metric functions have the limiting forms 
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2y/E' 



(72) 



The nature of the numerical results can be seen in Figure 7 below: 

One can also attempt to solve analytically for the kink solutions. (Note that, up to 
scaling, there is a unique such solution in each dimension D = 4n + 4.) Let us again 
consider the eight-dimensional case, and change variables so that (|7ip becomes 



dsl 



,,2p 



dp 2 

Hp) 



+ e 2p f(p) (a, - A') 2 + e 2p dnj 



(73) 
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The specific choice of coordinate gauge in (|73p is one that often allows one to obtain explicit 
solutions in terms of rational functions of e p . 

In the present case, we find that the metric is Ricci-flat if 

(1 + 8/-2/ 2 )/ 



~ (8/2 + 12// + /2) 

and the function f(p) satisfies the equation 

2/(2/ 2 - 8/ - 1)/ + 3(/ - 2)/ 3 + 2(21/ 2 - 34/ + 2)/ 2 
+ 12/(17/ 2 - 26/ + 2)/ + 56/ 2 (/ - 1)(5/ - 1) = , 



(74) 



(75) 



where a dot denotes a derivative with respect to p. This equation can be reduced to a single 
first-order differential equation as follows. We define f(p) = x and f(p) = y(x), which 
implies 



and ([75]) becomes 

3(x - 2)y 3 + 2(2x 2 



y 



y( z ) 



34x + 2)y 2 + 12x(17x 2 - 26x + 2)y + 56x 2 (x - l)(5x - 1) 



2x(l + 8x - 2x 2 )y 



(77) 



where the prime denotes a derivative with respect to x. Note that x is the squashing 
parameter, ranging from x = 1 near the origin of the kink, where the 7-spheres are round, 
to x = 1/5 in the asymptotic region near infinity, where the 7-spheres approach the squashed 
Einstein metric. The function h appearing in the metric (1T3f) is given by 



h 



x(l + 8x - 2x 2 ) 
8x 2 + 12xy + y 1 



(78) 



The kink solution corresponds to y(x) in (|77p describing an arc, lying below the y axis, 
starting at (x,y) = (1/5,0) and ending at (x,y) = (1,0). The asymptotic forms at the two 
endpoints of the arc are 



y 
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(5X_1) "I05 (5X " 1) + 5145 (5X_1) + 



24 



.2(1 - X ) + — (1 - X ) 



640 
343* 



'1 -xf + - 



(79) 
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Figure 7: The left-hand plot shows the functions a (lower) and b (upper) for the kink solution 
written using the metric ansatz (|7ip . The right-hand plot shows a 2 /b 2 , which ranges from 
1 (the round S 7 ) at small r to 1/5 (the squashed Einstein S 7 ) as r goes to infinity. 
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respectively. In terms of the original variables, the kink solution runs from 

/ ~ 1 , h ~ i (80) 

near the origin at p = — oo to 

1 9 

/ ~ - , h ~ — (81) 
^ 5 20 v ; 

as p approaches infinity. 

It should be remarked that there exists a known solution to the equations, corresponding 
to the complete metric of Spin(7) holonomy found in [19] I20j. This corresponds to 

/(p) = i(l-e- 10 ^ 3 ), or y(x) = ~(l-5x). (82) 

This solution has the same behaviour at large p (i.e. x — > 1/5) as we require for the 
kink solution, but it is very different at short distance (corresponding to p = in this 
parameterisation), since it has an S 4 bolt. In fact in the Spin(7) solution the variable x lies 
in the range < x < 1/5. 

By making further transformations, one can cast (|77p into a standard form for an Abel 
equation of the first kind. First, we make use of the known solution given in (182j) . and define 
a new dependent variable v(x), related to y(x) by 

_ / 2x 2 - 8x - 1 _ 3 n-1 
V[X) \x 2 (5x-l)v(x) 2(5x-l)/ ' 1 j 

After this change of variable, equation ([77]) becomes 

x 3 (5x - l)vv' = xV + 6x(4x 2 + 5x - l)v + 28(1 - x)(2x 2 - 8x - 1) . (84) 

The further change of variable to u(x) = (5x — l)~ 1 ^ 5 v(x) yields 

x 3 (5x - I) 7 / 5 uv! = 6x(5x - l) x / 5 (4x 2 + 5x - 2)u + 28(1 - x)(2x 2 - 8x - 1) . (85) 

It seems not to be possible to carry this further, since the change of independent variable 
required to put the equation into the canonical form udu/dz — u = q(z) is given by taking 

6(31x-2) _ 5 4/5 -1/5 p n 1.6. ± ) (m) 

X (5 X _ 1)1/5 ^ Z ° X 2r H 5 '5'5'5x^' y°°> 

yielding 

du _ 14(1 - x)(2x 2 - 8x - 1) 
U dz ~ U ~ 3x(5x - 1) 1 /5(4 X 2 + 5 X _ 1) " ^ ^ 



Since one cannot invert (|86p explicitly to obtain x as a function of u, it appears that no 
further progress towards an analytic solution can be made. 

Although attempts to solve (|77p completely by analytic means have not proved success- 
ful, one can use this first-order system to perform a phase-plane analysis. This can be done 
by writing the first-order equation (|77p in terms of an auxiliary parameter t, with 

^ = 2x(l + 8x-2x 2 )y, 
at 

j- = 3(x - 2)y 3 + 2(2x 2 - 34x + 2)y 2 + 12x(17x 2 - 26x + 2)y 

+56x 2 (x - l)(5x - 1) . (88) 
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The kink solution lies within the region 1/5 < x < lin the phase plane. In fact it 
corresponds to the unique flow that starts at the saddle at (x,y) = (1,0), and ends at 
the attractor at (x,y) = (1/5,0). It is also instructive to look at the exact Spin(7) solution 
([82]) . This starts at the saddle at (x,y) = (0,2/3), and flows (along a straight line) to the 
attractor at (x,y) = (1/5,0). The kink solution, and the extrapolation of the exact Spin(7) 
solution into the region 1/5 < x < 1, are shown in Figure 8 below. 




Figure 8: The numerical solutions of (177p corresponding, in the left-hand plot, to the exact 
solution of Spin(7) holonomy given in (|82p (extrapolated into the region 1/5 < x < 1 
depicted here), and in the right-hand diagram, to the kink solution. The latter corresponds 
to a trajectory from the saddle at x = 1, y = to the attractor at x = 1/5, y = 0. 

3.5.2 Stability of the kink soliton 

In the case of the kink the eigenvalue problem (I55j) for the spectrum of small perturbations 
around a static solution reduces to the single equation 

- ^e 5 ° (e*VV A )' + Kv x = X 2 v x , (89) 

where 

K = -L^-72e- 3 ' Bo -16e 4Bo +200e- 200Bo ) (90) 
+ 1 ( 2 88e- 3Bo - 48e 4B ° - 240 e - 200i?o ) B' 
+ I ( _ 288e~ 3B ° - 36e 4B ° + 72 e - 200B °) B' Q 2 . 

We have shown above that the profile function of the kink is monotone, hence the zero 
mode corresponding to the scaling freedom, vq = rB' (r), has no zeros which implies by the 
standard Sturm-Liouville theory there are no negative eigenvalues. Thus, the kink solution 
is linearly stable within our ansatz. 

The behavior of the kink under non-infinitesimal perturbations was studied numerically 
by the methods described in Section 3.4.2. As in the case of other static solutions we found 
that for small perturbations the kink is asymptotically stable, while for large perturbations 
it collapses to a black hole. The numerical evidence for these properties is shown in Figs. 9 
and 10. 
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Figure 9: Asymptotic stability of the kink solution. For small perturbations of the kink we plot a 
series of snapshots of (a) the function B(t,r) and (b) the function A(t,r). The dashed line shows 
the unperturbed kink solution. 



4 Conclusions 

In this paper, we have studied gravitational solitons and kinks in higher dimensions. Our 
focus has been the study of their stability, principally in the case of solitons in nine dimen- 
sions. We have considered various possibilities for the spatial metric, including examples 
such as the higher-dimensional Taub-NUT metrics [131 114j . which are not supersymmetric, 
and also the example of the As 8-metric of Spin (7) holonomy [151116] . which is supersymmet- 
ric. All the solitons we consider are trivial topologically (i.e. M. n topology), but non-trivial 
geometrically. We studied the question of stability first at the linearised level, using analyic 
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Figure 10: Instability of the kink for large perturbations. For large perturbations of the kink we plot 
a series of snapshots of (a) the function B(t,r) and (b) the function A(t,r). During the evolution 
A(t, r) drops to zero at r = rji ~ 1.754 which signals the formation of an apparent horizon there. 
Outside the horizon the solution relaxes to a static black hole. 



methods, and found that in all cases the solitons are linearly stable. Numerical analysis 
indicates that this stability persists for non-infinitesimal perturbations also, provided they 
are small enough in magnitude. The numerical analysis also shows that for sufficiently large 
perturbations the solitons are all unstable to black-hole formation. This is not altogether 
surprising, since even flat Minkowski spacetime is unstable to sufficiently large perturba- 
tions, which can lead to the formation of a black hole. These instabilities provide a salutory 
reminder of the fact that in gravitational theories supersymmetry is not a guarantor of 
stability beyond the linearised level. 

We also studied in detail a kink metric in eight dimensions, whose existence was first 
encountered in [18]. It is a non-trivial cohomogeneity-1 metric on R 8 , in the which the 
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level sets are homogeneous 7-spheres viewed as S 3 bundles over S" 4 . At small distances the 
level surfaces approach the round 7-sphere and the metric is of the form near the origin of 
Euclidean 8-space. At large distances the level surfaces approach the squashed S 7 Einstein 
metric in the family of S 3 bundles over S . We showed that the Ricci-flatness conditions for 
the kink metric can be reduced to a first-order Abel equation of the first kind, but it appears 
not to be possible to obtain an explicit solution analytically. Our discussion includes a proof 
of the existence of the Ricci-flat metric. 
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A Nine-Dimensional Schwarzschild 

This section contains a brief summary of the analysis of the stability of the nine-dimensional 
Schwarzschild solution within the framework of the deformations considered in this paper. 
If we expand around the nine-dimensional Schwarzschild background, for which 

A = l-^, 8o = 0, ^ 0) = 0, 4° ) = 0, (91) 

we find that working to first order in fluctuations we can keep A = Aq and 5 = 0. For 
the dynamical fields, we shall now use c/>i and <p 2 to denote the linearised fluctuations. 
Introducing the "tortoise coordinate" x via dx/dr = A^ 1 , and defining 



J x = r' 7 / 2 ui, 02 = r- 7/2 u 2 , (92) 



we find that these satisfy 



ili--^ + V Ul =0, U2--^- + Vu 2 = 0, (93) 
where the potential V (which is the same for both u\ and u 2 ) is given by 

"?)(")• <«> 

Note that this is the same potential as was encountered in |17j in the analysis of the nine- 
dimensional perturbations of Schwarzschild with a single dynamical variable (corresponding 
to u 2 = here.) 



B Time Evolution of Kerr-AdS Black Holes 

In this appendix, we use the techniques of this paper to analyse the stability of a particular 
class of five-dimensional rotating black holes. The general methods extend extend to any 
black hole having all angular momenta equal, but in this appendix we restrict attention 
to the (4 + l)-dimensional case. This means that we can again consider an ansatz with 
SU(2) x U(l) isometry on the constant-radius spatial sections. 
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The five-dimensional Kerr-AdS solution [21] with equal angular momenta is given by 
dsl = - . ° d t^ + d 4r + \P 2 f 1 + (-3 + 2^ dt) 2 + \p 2 (a 2 + 4), (95) 



2Ma2\ G Ar V p 4 



+ ^ 

where 

/ 2 2MH 2Mo 2 
G = 1-Ap 2 — + 

The metric satisfies fi^Af = 4A<7mn- Thus for the asymptotically AdS case we require 
A < 0. 

The new feature in our analysis is the inclusion of the Kaluza-Klein vector in the di- 
mensional reduction. The electric charge associated with this field is proportional to the 
angular momentum of the black hole. Thus we make the following reduction ansatz 



ds 2 = e 2a ^ds 2 + e 2 ^ 



^^a z + A) 2 + e~ 2 ^{a 2 + a z 2 ) . (97) 



For the time being, we consider the base metric ds 2 to have dimension n. Later, we shall 
specialise to the case of immediate interest, namely n = 2. We choose the natural vielbein 
basis 



e° = e^ +2 ^(^+A), 



e a = e a , 



a = 1,2. (98) 



After some calculation, we arrive at the following non-vanishing Ricci-tensor components 
(in the vielbein basis): 

Roo = e- 2 ^{-(3/3+(n-2)a)[/3(Vv9) 2 +2 7 V^-V^-/3n^-2 7 n(/>} 

R a /3 = e- 2Q ^{(6a/3-3/3 2 + (n-2)a 2 )V Q ^V /3V 9-67 2 V Q( /)V^-(3/3+(n-2)a)V Q V^} 
-a e - 2 ^{(3/3+(n-2)a)(V^) 2 +D^}^ 
e~ 2aip |(3/3+(n-2)Q)[ 7 Vv?-V^-/?(V^) 2 ]+7D ( /.-/3n(p} 



R, 



_l e -2/3^+2 79 i ( e 6 7 ^_ 2 ) 



R 0a = ^-P^-^y^iW-^+A^ ^A + iina+S^e^-^+^T^Vpip, (99) 
where T = dA. 
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Specialising to the case n = 2, these expressions give 

R 00 = e ~ 2a v> {-3(3[(3{Vtp) 2 + 2~fVip-V(f)]-0\3ip-2~(n0} 

+ l e (2/3-4a)^+4 7 l e -2^+8 7 ^ (10Q) 

R af 3 = e- 2 ^{3/3(2a-/3)V Q <^V^v9-6 7 2 V cl 0V^-3/3V Q V / 3^} 

-ae" 2 ^ {3/3(V 95 ) 2 +n^}7 ?a/3 -i e ^- 4 ^ +4 ^^ 7 ^+e- 2 ^ J R a/3 , (101) 
e -2^| 3 ^[ 7 Vy,.V0-^(V^) 2 ]+7n^-/3n^}-ie- 2 ^ +2 ^ (e 6 ^-2)] S a (102) 

1 e -(a+4/3) V -2 7 ^ ^(5/3-2 a ) V +4 7( /, ^ _ (1Q3) 



Following the analysis in [Sj, we use the volume of the S 3 to parameterise the radial 
direction, and write 

d f = _ Ae -2S df+A- 1 dr 2 + \r 2 [e 2B (al+al)+e- iB (a 3 +A) 2 }. (104) 



Thus, comparing with (|97[) . we have 

e^ = \r, 70 = -B, ds 2 = e~ 2a ^ (- Ae~ 2S dt 2 + A~ l dr 2 ) . (105) 
Substituting into the Einstein equation 

Rmn = 4A gMN ) (106) 

we note from f|103[) that 

V = c e „e( 2a - 5 ^, (107) 

where c is a constant and e^ u is the Levi-Civita tensor in the metric ds 2 . This expression 
can be substituted into the remaining Einstein equations, leading to the momentum and 
Hamiltonian constraints 

A! = -^ + L{Se- 2B -2e^ B )-2r{e 2& A~ 1 B 2 +AB' 2 )- l ^e iB -A\r, 
r 3r 3r ' 

A = —ArABB' , (108) 
the slicing condition 

5' = -2r(e 2S A~ 2 B 2 +B' 2 ), (109) 

and the wave equation 

,<5 A-l„3 dY ( „-S A „3 d/V i 4 -(5„/< -2B -8B\ , 128c 2 j 4B 



eM-^Sj -[e-°Ai#B') +^e~ 6 r ( e -^- e -«*) + rp_ e -« e *« = . (110) 

It is straightforward to verify the self consistency of the equations. Namely, that the 
vanishing of the dot of A 1 minus the prime of A in (|108p yields, after using (|109|) and (|108p . 
the wave equation (|110|) for i?. 

Defining j4 = 1 — m(r,t)/r 2 — Xr 2 , the Hamiltonian constraint in (|108p becomes 

m ' = 2r 3 (e 25 A- 1 B 2 +AB' 2 ) + lr(3+e- 8B -ie- 2B ) + ^^e AB . (Ill) 
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This is manifestly positive. 

Note that the constant c is related to the angular momentum. Lowering the index on 
the Killing vector d/dtp gives the 1-form 

K = f&v+W (a 3 +A) . (112) 

The angular momentum is given by the Komar integral J = l/(167r) J s3 *dK, where * is 
the Hodge dual in the five-dimensional metric ds 2 , and we have 

*dK = e (5/?-2«) V +4 7 4> *jFA(7i Aa 2 Aa 3 + - ■ ■ , (113) 

where *J- is the 0-form Hodge dual of the field strength T in the 2-metric ds 2 . Thus the 
angular momentum can be seen to be nothing but the 2-dimensional electric charge of the 
dimensionally-reduced solution. From (|107[) . the Komar integral therefore gives the angular 
momentum 

J = -77- CT1A0-2A0-3 = vrc. (114) 

l07T J S 3 

Comparing the radial coordinate r used in (|lU4p with the radial coordinate p used in 
([95]) . we see that 

r = K 1+ ~^~) ■ ( 115 ) 

This can be used to rewrite the previous equations using p instead of r as the radial variable. 
We make the expansion 

A = A (1+A), B = B +B, 5 = 5 +5, (116) 

where Aq, Bq and 5q denote the background expressions in the Kerr-AdS metric, which can 
be read off from (|95p and (|96|) . and we work to linear order in the tilded quantities. 
From the slicing equation (|109j) we obtain 

(3p 4 + 2Ma 2 )^ + 16Ma 2 ^ = 0, (117) 
op op 

whilst from the Hamiltonian constraint in (|108j) we find 

(3p 4 + 2Ma 2 )A+16Ma 2 B = 0. (118) 

(Taking a constant of integration to be zero.) These two equations can be used to solve 
for the perturbations 5 and A in terms of the dynamical variable The momentum 
constraint in (|108|) implies 

(3p 4 +2Ma 2 )A + 16Ma 2 B = , (119) 

which is consistent with (|118p . 

The wave equation (jllOp then gives 

YlMa? 1 16M 2 a 2 (3+2Aa 2 ) 4M 2 a 4 



p 3 / , 2Ma 2 \~ . 3 , 8p(3+^+ 3p , 



1+^ B-(^GB')'+ r - ^ '-^B = 0, (120) 



7 Note that in the case of linearisation around the Schwarzschild solution, discussed in [8], one can take 
5 = and A = 0, so that only the perturbation B of the dynamical variable B need be considered in that 
case. 
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where a prime here means d/dp, and, from <\96\ 



, 9 2ME 2Ma 2 , . 

G^l-Xp 2 — + ^^- 121 

We can cast the wave equation into a Schrodinger form by introducing the "tortoise 
coordinate" x defined by 

cix=(l + — G dp, (122) 

and introducing a new dynamical variable x(a?,t), defined by 

2Ma 2 \-i/4 
P 

Equation (|120p then takes the form 



S = p-3/2^ 1 + ^_j x . (123 ) 



*-0 + y ^ = ' ( 124 ) 



where the potential 1/ is given by 
V 



G 

v 



, o 160M 3 a 4 (l+2Aa 2 )p 2 16M 2 a 2 (7p 2 + 15Aa 2 p 2 +a 2 
35+15A/9 +- 



(p 4 + 2Ma 2 ) 3 (p 4 +2Ma 2 ) 2 
2M(87p 2 + 58Aa 2 p 2 - 92a 2 ) 256M 2 a 2 (3p 2 + 2Aa 2 p 2 - 2a 2 ) 



+ 



p 4 +2Ma 2 (3p 2 +2Ma 2 ) 2 
64M(9p 2 + 6Aa 2 /j 2 -8a 2 ) 



3p 2 + 2Ma 2 



(125) 



The structure of the potential V can easily be studied in the special case when A = 0, 
so that the background metric is the Ricci-flat a = b Myers- Perry solution |22j . It can be 
seen from (|12ip that in order to have an horizon, i.e. for the function G to have a zero for 
real p, it must be that a 2 < \M. If, therefore, we write M = 2a 2 (l+s 2 ), where s is a real 
constant, then the outer horizon occurs at 



p\ = 2a 2 (l+s 2 + s\/l + s 2 ) • (126) 

Writing p 2 = p^+y 2 , and substituting this and M = 2a 2 (l+s 2 ) into the expression (|125|) for 
the potential, one finds that V is manifestly non-negative everywhere outside the horizon, 
and it tends to zero on the horizon and at infinity. 
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